Avoided Quantum Criticality near Optimally Doped High Temperature 

Superconductors 
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We study the crossover from the underdoped to the overdoped regime in the t-J model. The un- 
derdoped regime is dominated by the superexchange interaction, locking the spins into singlets which 
weakly perturbe coherent charge carriers. In the overdoped, large carrier concentration regime, the 
Kondo effect dominates resulting in spin-charge composite quasiparticles which are also coherent. 
Separating these two Fermi liquid regimes, there is a critical doping where superexchange and Kondo 
interaction balance each other, bringing the system close to a local quantum critical point near the 
point of maximal superconducting transition temperature. At this point, particle hole symmetry 
is dynamically restored and physical quantities such as the optical conductivity, exhibit power law 
behaviour at intermediate frequencies as observed experimentally. Quantum criticality is avoided 
by the onset of superconductivity. 

PACS numbers: 71.27. +a,71.30.+h 



The high temperature superconductors have a normal 
state which is not well described by the standard theory 
of metals pj . Several anomalies are well established ex- 
perimentally: a) The normal state close to the doping 
level with the highest transition temperature is highly 
incoherent and does not fit the standard Fermi liquid 
picture, b) at small doping the normal state has a pseu- 
dogap, c) the optimally doped regime exhibits anomalous 
power laws in transport quantities and in the optical con- 
ductivity in an intermediate frequency range. This obser- 
vation has lead to the suggestion that the superconduct- 
ing dome covers a critical point , d) The superconduct- 
ing state not only exhibits coherent Bogolubov quasipar- 
ticles and charge collective modes, but also a very sharp 
spin mode or " 40 meV resonance" . The role of this mode 
as well as of the putative quantum critical fluctuations 
in a fundamental theory of cuprate supercondcutivity are 
important issues. 

In this letter we present a coherent account of these 
observations based on a Cluster Dynamical Mean Field 
Theory analysis of the t-J model. This approach allows 
us to access the regime of intermediate temperatures in 
the doping regime between underdoped to overdoped sys- 
tem. In addition, the approach is formulated in terms of 
pseudoparticles, representing collective excitations of the 
system in an enlarged Hilbcrt space. This formulation 
provides a clear interpretation of the numerical results, 
providing a transparent physical interpretation. 

We start from the the t-J model which was proposed 
by P.W. Anderson as a minimal model to describe the 
cuprate superconductors 
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It contains two terms; one describes the kinetic energy 
which delocalizes the holes introduced by doping, and a 
spin spin interaction. A constraint of no double occu- 



pancy must be enforced. 

We used the extended Dynamical Cluster Approxima- 
tion Q combined with the slave particle approach £| 
which have been very useful in the theory of strongly 
correlated electrons. The cluster approach maps the full 
many body problem, into a a set of local degrees of free- 
dom (in this case a two by two plaquette) which are 
treated exactly and a bath which is treated self consis- 
tently. 

The coupling of the cluster to the medium, which sim- 
ulates the rest of the lattice, causes the cluster eigen- 
states to decay in time therefore their spectral functions 
carry nontrivial frequency dependence and important in- 
formation about various physical processes such as the 
RKKY interactions, the Kondo effect and d-vawe super- 
conductivity. To study the dynamics of the cluster states, 
we diagonalized the cluster, i.e., i? c ; us t er |m) = E m \m), 
and assigned to each cluster eigenstate a pseudoparti- 
cle m) = aJjjO). The original problem can be exactly 
expressed in terms of pseudoparticles a m with the only 
non-quadratic term of the converted Hamiltonian being 
the hybridization between the cluster and the medium. 
This cubic hybridization term is taken into account by 
the self-consistent resummation of the non-crossing dia- 
grams in the spirit of the well known Non-Crossing ap- 
proximation to the Anderson impurity model. Since the 
Hilbert space of the new base of pseudoparticles accesses 
some non-physical states, one needs to project them out 
with the constraint expressing the completeness of the 
cluster eigenbasis. First we calculate the pseudopar- 
ticle self energies and the pseudoparticle Greens func- 
tions. Finally, physical observables like spectral function 
or susceptibilities can be calculated by a weighted sum of 
the convolution of the pseudoparticle spectral functions. 
The following matrix elements (F^) nm — (n|c^-|m) and 
(Sq)™ = (n|SQ|m) give the bare weight of the electron 
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FIG. 1: The cluster (tt,0) self-energy at zero frequency as 
a function of temperature for few doping levels. The inset 
shows the estimation of the coherent scale in the normal state 
of the t-J model (black dots) and transition temperature to 
superconducting state (red dots). 

spectra and the spin response in terms of pseudoparticles 

GkM = - ( FK )rn>n>C m >„> nm (iu)(F K % m (2) 

nmn'm' 

Xq {i^) — ^ ' (SQ) m i n iC m i n i nm (iuj)(S^_Q) nm (3) 

nmn'm' 

with 

C m >n' nm (i<j) = T^G n > n (ie)G mm > (it - iw) (4) 

ie 

Here Gk, G nn i and xq are electron Green's function, 
pseudoparticle Green's function and spin susceptibility, 
respectively. This is the central equation of the approach, 
relating observables to the pseudoparticle spectral func- 
tions plotted in Fig.Q] More details of the method can be 
found in Q . It has been shown that the cluster approach 
successfully describes many properties of the high tem- 
perature superconductors y. Li li| In this letter we show 
that it also unravels the origin of the apparent criticality 
observed near critical doping, and its elimination below 
the superconducting state. 

Physically pseudoparticles represent coarsed grained 
versions of the important many body excitations includ- 
ing fermionic quasiparticles and bosonic collective modes. 
They have quantum numbers describing their spin, num- 
ber of particles, (which divided by the cluster size, give 
the density), and a coarsed grained momentum. In mo- 
mentum space, the course graining for the smallest clus- 
ter divides the Brillouin zone into four different regions, 
namely an inert patch around (0, 0) point, the second 
patch centered around (jr, tt) being highly non-Fermi liq- 
uid like at small doping and finally (jr, 0) and (0, n) re- 
gions which contain the Fermi surface of the model in 
the underdoped and slightly overdoped regime and are 
therefore most relevant for thermodynamical and trans- 
port properties. 
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FIG. 2: a) The absolute value of complex optical conduc- 
tivity <t(u) is proportional to a;" 2 / 3 in the intermediate fre- 
quency region for optimally doped system, b) The absolute 
value of the (tt, 7r) spin sussceptibility is proportional to u~ 2 . 
c) The angle, calculated as arctan(Im(<j)/Re(o-)) is approxi- 
mately 7r/3 since <t(lu) ~ (— iio)" 2 ^ 3 . In the legend, N stand 
for normal state and S for superconducting state. 

The first indication for an underlying criticality near 
optimal doping comes from the evaluation of the electron 
scattering rate obtained from the imaginary part of the 
electron self-energy as a function of temperature for few 
different doping levels. As described in Fig. [T\ both at 
large and small doping the scattering rate is small as ex- 
pected for a Fermi liquid. Remarkably it becomes very 
large in the region near optimal doping when the critical 
temperature is maximal. The transition to the supercon- 
ducting state severely reduces the scattering rate elim- 
inating the traces of the underlying critical behaviour, 
hence the name avoided quantum criticality. A coher- 
ence scale, estimated from the scattering rate, is plotted 
in the inset of Fig. [T\ and shows it tends to vanish close 
to the point of maximal superconducting transition tem- 
perature. 

Additional evidence for the quantum criticality is the 
emergence of power-law behaviours of the response func- 
tions. As shown in Fig.El the optical conductivity has an 
approximate power law with an exponent 2/3. The same 
powerlaw is realized in the one-electron self-energy while 
the spin susceptibility is proportional to Xw,n °c lu~ 2 
in the same frequency range. Experimentally, it was 
found that the optical conductivity is proportional to the 
(— iuj)~ - 65 in the intermediate frequency regime 0. 

Quantum criticality is avoided when the electrons con- 
dense forming d-wave pairs. The electron scattering rate 
is dramatically reduced (see Fig. and a V-shaped gap 
opens in the local one-electron density of states. The 
particle-hole response at (tt, tt) is severely reduced for 
frequencies below the superconducting gap and a sharp 
resonance appears in the gap. A sharp resonance peak in 
the spin susceptibility was observed in bilayer cuprates 
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FIG. 3: The dynamical spin susceptibility at q = (it, tt) for 
few different doping levels and three different temperatures: 
superconducting state, normal state at the transition tem- 
perature and at room temperature. The pronounced peak 
is formed in SC state at Q.lQt ~ 48meV and a broad peak 
in normal state is around 100 — 140 meV. Susceptibility at 
normal temperture is much smaller and the peak moves to 
higher frequencies. The resonance is strongest at the opti- 
mally doped system. It disappears quickly in the overdoped 
site and somewhat more slowly in the underdoped side. 
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FIG. 4: left: Pseudoparticle spectral functions for the three 
most important pseudoparticles: ground states for N=4, N=3 
and N=2 sectors. Right-top: Sketch of pseudoparticle tresh- 
old energies which can be interpreted as the effective many- 
body levels in normal and superconducting state. Right- 
bottom: Pseudoparticle occupancies versus doping for most 
important pseudoparticles. The full lines correspond to the 
normal state while the dashed lines correspond to the super- 
conducting state. 



YBCO and Bi2212 and in singe layer T12201. It scales 
with doping like 5T C and does not depend on tempera- 
ture. 

Our method gives a pronounced peak in the (n, 7t) 
spin response at frequency 0.16i in the optimally doped 
regime, as shown in Fig. EI The peak position is temper- 
ature independent and taking into account the maximal 
transition temperature in our model T c - max ~ 0.035< the 
resonance energy is 4.6T C in good agreement with exper- 
iments. 

In addition we see a broader peak around 0.35 — 0A5t 
extending to very high frequencies of order of t « 
300 meV which also gains some weight upon condensa- 
tion and represents an important contribution to the ex- 
change energy difference between superconducting and 
normal state. This exchange energy gives far the largest 
contribution to the condensation energy 

Our method averages the momentum dependence over 
1/4 of the Brillouin zone centered at (tt, 7r) therefore it is 
reasonable to compare our results with the q integrated 
susceptibility from Ref. where both the 35 meV reso- 
nant peak as well as broader peak around 75 meV ex- 
tending up to 220 meV was observed. 

Pseudoparticle interpretation: These results can be un- 
derstood in terms of the evolution of the pseudoparticle 
spectral functions (spectral functions of the eigenstates 
of the plaquette) depicted in Fig. 0] Remarkably, out of 
the large number of pseudoparticles which are involved 
in our calculation (3 4 ) only four distinct pseudoparti- 
cles are very important, having the largest weight in 
the ground state and containing more than 95% of the 



spectral weight of the low energy correlation functions. 
Fig. 0k shows the evolution of the three most important 
pseudoparticle spectral functions from the underdoped to 
the overdoped regime. 

Underdoped Regime: At small doping, the singlet 
state with one particle per site and zero momentum 
(N=4,S=0,K=0) (half filled singlet) dominates since it 
has far the lowest treshold energy, is very strongly peaked 
at the treshold energy, and has the largest occupancy as 
shown in Fig.QJ:,. It describes a system locked in a short 
range singlet state as a consequence of the strong su- 
perexchange interaction. The electron spectral function 
descrbies the process of addition and removal of an elec- 
tron from the system at frequency to. It is dominated by 
convolution of two pseudoparticles with different cluster 
occupation N and N+l, or N-l and with integral in con- 
volution taken between zero and u>. The necessary con- 
dition for a peak of the one-particle spectral function at 
the Fermi level is that at least two pseudoparticle spec- 
tral functions share common threshold and are strongly 
peaked at the same threshold. Since the thresholds of 
the other pseudoparticles are significantly shifted with 
reference to the half filled singlet, a pseudogap results in 
the one particle spectra in the underdoped regime. This 
gap in threshold energies severely limits the possible de- 
cay processes of the electron resulting in a low electronic 
scattering rate. 

The local picture of the overdoped regime: On the over- 
doped side, all three important pseudoparticles (half- 
filled singlet, doublet with one hole per plaquette and 
singlet with two holes per plaquette) have a power law 
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divergence at the same threshold frequency at zero tem- 
perature (Fig. |2k and Fig. ^t) which is a standard sig- 
nature of the Kondo effect. Hence, the one particle spec- 
tral function develops the Kondo-Suhl resonance at the 
Fermi level since the convolution between the doublet and 
half-filled singlet (or N=2 singlet) state is large at low 
frequency. The one-particle spectral function is peaked 
slightly above the Fermi level. The underdoped regime 
is distinctively different since only the half-filled singlet 
state dominates while the doublet has very little spectral 
weight in the region of the singlet peak. 

The Transition Region, Normal State: In the optimal 
doped regime, the Kondo effect and the superexchangc 
compete giving rise to a regime with very large scattering 
rate and consequently a small coherence scale. Criticality 
is the result of the merging of the convergence of the 
thresholds of the pseudoparticles (See Fig 4) when going 
from the underdoped to the overdoped regime. Local 
criticality has been found in models of heavy fermion 
systems Q and frustrated magnets . 

Surprisingly this evolution is accompanied by a 
restoration of particle hole symmetry in the optimally 
doped regime. As shown in Fig. the threshold of the 
N=2 cluster ground state and N=3 cluster ground state 
(doublet) merge first resulting in a Kondo-likc contribu- 
tion to the electron spectral function. This contribution 
is peaked above the Fermi level as usual for the sim- 
ple Fermi-liquid state in less than half-filled one band 
model. The half-filled singlet however remains the lowest 
state in energy and still gives a significant contribution to 
the electron spectral function. The later contribution is 
peaked below the Fermi level and keeps a pseudogap-like 
shape. The result of the two contributions to the elec- 
tron spectral function is a restoration of the particle- hole 
symmetry in the density of states both in normal and 
superconducting state at optimal doping. The particle- 
hole symmetry of the density of states is yet another hint 
of the proximity to the quantum critical point since it is 
well known that only this symmetry allows non trivial 
criticality in the two impurity Kondo model 

Transition into the superconducting state: The degen- 
eracy responsible for the strongly incoherent metal with 
large scattering rate at the Fermi level is lifted by the 
superconductivity avoiding the quantum critical point. 
Fig. shows that both important singlet pseudoparti- 
cles (for N=4 and N=2) develop a very sharp peak at 
the same treshold frequency and, at the same time, their 
occupancy increases (see Fig.Hfc) upon condensation, in- 
dicating that electrons are locked into singlets with zero 
momentum. A gap opens between the singlets and dou- 
blets which gives the gap in the one-particle density of 
states. Because of this gap in the pseudoparticle thresh- 
olds, the large imaginary part of the electron self-energy 
does not persist in the superconducting state (see also 
Fig. ^| . Since the density of states is composed of two 
almost equally important contributions, i.e., the convo- 



lution of the doublet with both singlets (N=4 and N=2), 
the superconducting gap is almost particle hole symmet- 
ric in the optimally doped regime with half-width of the 
order of 0.12t. When the doping value is changed from its 
critical value, the asymmetry in the superconducting den- 
sity of states appears. The magnitude of the assymetry is 
the same as the asymmetry of the corresponding normal 
state spectra and comes from the fact that the occupancy 
and therefore importance of the N=4 half-filled singlet 
exceeds the importance of the N=2 singlet (see Fig. Efc) . 
The spin susceptibility comes almost entirely from the 
convolution of the half-filled singlet with the half-filled 
triplet (N — 4, S — 1,K = (tt,7t)). The later develops a 
peak at an energy 0.16i upon condensation which causes 
the resonance in the spin susceptibility. 

In conclusion we showed that at optimal doping, in the 
t-J model, there is a collapse of coherence scale, accompa- 
nied by a a large scattering rate, restoration of particle- 
hole symmetry, which result in approximate power laws 
in physical quantities, such as the optical conductivity 
(a oc cj~ 2 / 3 ) and spin susceptibility (xu-tt cx uj~ 2 ). Quan- 
tum criticality is avoided when electrons condense into 
a supercondacting state which results in a dramatic re- 
duction of scattering rate, the formation of an energy 
gap which in turns leads to a resonance in the spin 
susceptibility. A transparent physical interpretation of 
the numerical results was given in terms of pseudopar- 
ticles which represent the important many body exci- 
tations, fermionic quasiparticles and bosonic collective 
modes. This real space picture complements the momen- 
tum space picture which emerges in other studies, where 
substantial deformations of the Fermi surface which re- 
sult in the formation of Fermi arcs at comparable values 
of the doping 0, |^. Finding analytic connection be- 
tween single particle properties in k space and the collec- 
tive excitations presented in this study is an interesting 
direction for future research. 
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